Complete Condensation in a Zero Range Process on Scale-Free Networks 
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We study a zero range process on scale-free networks in order to investigate how network structure influences 
particle dynamics. The zero range process is defined with the rate p{n) = n^ at which particles hop out of nodes 
with n particles. We show analytically that a complete condensation occurs when 5 < 5c = 1/ (y— 1) where y is 
the degree distribution exponent of the underlying networks. In the complete condensation, those nodes whose 
degree is higher than a threshold are occupied by macroscopic numbers of particles, while the other nodes 
are occupied by negligible numbers of particles. We also show numerically that the relaxation time follows a 
power-law scaling x ~ U with the network size L and a dynamic exponent z in the condensed phase. 
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Network systems have been attracting much interest in re- 
cent years 1 1, 2, 3, 4]. A network is a simplified description 
for complex systems, such as the Internet, with a set of nodes 
interconnected with links. Various networks have been intro- 
duced and studied in wide areas of science, many of which 
share the common feature of the power-law degree distribu- 
tion Pdeg. W ~ A:-^ f?) . A degree of a node denotes the num- 
ber of links attached to it. Such a network is called a scale- 
free (SF) network. 

The SF network has the broad degree distribution, which 
leads to many intriguing properties of dynamical systems de- 
fined on it. One of those is a centralization. For instance, a 
random walker on a network tends to visit more frequently 
those nodes with higher degree |,5]. Higher degree nodes usu- 
ally have large values of the load or betweenness centrality, 
which is a measure of a traffic along shortest pathways |6,|7|]. 
Whether under stochastic or deterministic dynamics, particles 
tend to concentrate on high degree nodes. Consequently, one 
expects that a particle interaction would play an important 
role in dynamic processes. However, only a little attention 
has been paid to interacting dynamical systems on SF net- 
works f?!. 

In this work we study a zero range process (ZRP) on 
SF networks. The ZRP is a model for condensation in inter- 
acting stochastic particle systems. It has been widely studied 
since it has a lot of applications such as traffic flows and its 
stationary state properties are solvable. In this study, we in- 
vestigate the interplay between particle interaction and under- 
lying network structure. 

Consider a ZRP for particles on a network of L nodes. 
We mainly focus on the limiting case where N,L ^ °o with 
fixed density p = N/L. Each node / = l,...,L may be oc- 
cupied by any number of particles. We denote the occupa- 
tion number at each node / with n, and a microscopic state 
with n = (ni ,n2, ■ ■ ■ jWl)- In dynamics, only a single parti- 
cle is allowed to jump out of each node during a unit time 
interval , and the jumping rate pi out of a node ; de- 

pends only on the occupation number at the source node. For 
simplicity, we assume that the jumping rate is given by the 
same function of the occupation number for all nodes, that 
is. Pi — p{n). A jumping particle out of / then hops to a target 



node j according to a prescribed transition probability denoted 

by Tj^i. 

The model incorporates particle interactions with the hop- 
ping rate p{n). When p{n) °^ n, particles are noninteracting 
and move independently. A repulsive interaction is imposed 
when p{n) increases faster than n. An attractive interaction 
is imposed when p{n) increases sublinearly [weak attraction] 
or deceases [strong attraction] with n. The ZRP may also 
model a transport process, where e.g., data packets are de- 
livered from one node to another. In that context, p{n) corre- 
sponds to the amount of packets that can be handled per unit 
time by each node with n packets, that is, the jumping rate 
function may represent a transport capacity of each node. 

The ZRP is solvable as far as stationary state properties are 
concerned. We give a brief review in the following. For de- 
tailed review, we refer readers to Refs. [9, Let fstat.(n) 
be the stationary state ^probability distribution. Interestingly, it 
has the product form 



(n) 
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where /,(«) = flLi [^'/f (01 for n > and /;(0) = 1 and 
Z{L,N) is a normalization constant. Here, {co,} is the sta- 
tionary state probability distribution of a single particle or 
a random walker moving on the network with the transition 
probability The mean occupation numberm,- = («,) — 

Ln «i-Pstat. (n) is given by m,- = z(3/3z)lnF,(z), where 



n=0 



(1) 



The series is defined within the interval |z| < Zc, where Zc is the 
radius of convergence. The fugacity z should be determined 
self-consistently from p = (Lf=i mi)/L. 

In the periodic lattices with CO,- = const., the condensation 
occurs when p{n) decays faster than pc{n) = 1 +2/n and 
p > Pt with a critical density p^ In the condensed state 
a node is occupied by macroscopic numbers of particles. It 
will turn out that the condensation occurs in a different way in 
SF networks. 
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Now we take a SF network of L nodes with the power- 
law degree distribution Pdeg.(^o l^k < ki„ax) ^ k^"^ ■ The de- 
gree distribution exponent is y and is a constant. We only 
consider the network with y > 2 so that the mean degree 
A: = / dkkPi^g {k) is finite. The node with the maximum de- 
gree kmux is called the hub. For a given Pi^g ik), the max- 
imum degree scales as k^a^ ^ which results from 

£„,fdeg.(feMfc = o(i/i). ' 

We adopt the following jumping rate function 

p{n)=n\ (2) 

where the parameter 5 controls the interaction strength or 
the transport capacity. A particle jumping out of a node / 
is allowed to hop to one of the neighboring nodes of / se- 
lected randomly. That is, the transition probability is given by 
Tj^i = l/ki if / and j are connected and 7}^, = otherwise, 
where fc, denotes the degree of 

When there is only a single particle, the ZRP reduces to 
random walks with the transition probability The ran- 

dom walk problem is studied in Ref. |5], where it is shown 
that the stationary state distribution is given by CO, — ki/E 
with E = 'Ljkj. One then obtains /(n) = n"=i = 
(ki/E)" /{n\)^. Hence the function in Eq. ((T) takes the 
form Fi{z) ^ !F^{x^ zhjE) with 



and the mean occupation number is given by 



ax 



(3) 



(4) 



x=xj=zki /E 

The parameters, = zkj/E will be called the effective fugacity. 
The higher the degree is, the larger the effective fugacity is. 
Note that the particle distribution is determined solely by the 
degree distribution. The mean occupation number is parame- 
terized with the fugacity z, whose value has to be fixed from 
p = (Lm,)A. 

First of all, we derive the mean occupation number dis- 
tribution in the simplest case with 5 = 0, where fo{x) = 
1/(1 — x) for X < 1 and m,- = Xi/{1 — x,). The actual occu- 
pation number is determined after solving the relation p = 
L^^ ^^^[ —Xi). In the L ^ oo limit, one might replace 
the summation with the integral. However, one has to be cau- 
tious since m as a function of x is singular at x = 1 . The hub 
has the largest effective fugacity, hence one may safely de- 
compose the particle density into two terms as p = p, + p„ for 
the hub and the other nodes, respectively. They are given by 



L 

kmax 
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dkPdegXk) 
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with k' = J^^'"" dkk'Pieg.{k)- Using the power-law degree dis- 
tribution, one obtains that k' ^ k^max ^/^ for large I. Hence 



the total density at all nodes but the hub is given by p„ ~ 
kmax ''110(1 — zkmax/E)\. One can eliminate the fugacity z 
from the relation for p.^, which yields that p„ ~ kmax ln(l + 
mhub)- Note that A;„K,x ~ i''^'^^'' and that ot/,,,^ < pL. It implies 
that p„ vanishes as p„ ~ L^' InL in the L ^ °o limit. Conse- 
quently, the whole fraction of particles is condensed into the 
hub with mfiuh /L~p. 

The analysis shows that a complete condensation occurs at 
5 = 0, where the whole fraction of particles is concentrated at 
the hub. Moreover, the condensation takes place at any finite 
value of p. With these two features the complete condensa- 
tion contrasts well with the usual condensation occurring in 
homogeneous networks. 

In the non-interacting case (5 = 1), it is trivial to show that 
m, = co,A^ = pki/k with the mean degree k. The hub also 
has the largest occupancy, but it scales sublinearly with L as 
minih ^ k,„ax ^ L'/'^^''. So the condensation does not occur 
at 5 = 1 . It will be shown that there exists a critical value of 5 
for the condensation. 

For the general cases with 5 > 0, we need to figure out the 
property of f^ix) in Eq. (|3}. The radius of convergence is in- 
finite and it is analytic for all values of x. Since a closed form 
expression is not available, we try to find an approximation. 
For large x, the series is dominated by the terms near n = no at 
which the summand x"/(n!)^ = e'''^"' is maximum. Using the 
Stirling formula and treating « as a continuous variable, we 
obtain no — x^^^ from g'{no) = 0. Then, the series is approxi- 
mated as f^{x) ~£~^oexp[§(no)+g"(no)(n-«o)^/2 + ---]. 
Replacing the summation with the integral, we obtain that 
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Note that the approximation coincides with the exact form 
(x) = g-* at 5 = 1 . We have checked numerically that the 

approximation is also good for 5^1 unless x is too close to 

zero. For small x, we can approximate the series with a few 

lowest order terms. 

Within the approximation, the mean occupation number is 

given by 



nii : 



X,- for X/ ^ 1 , 
x]^^ forx, > 1 . 



(6) 



The self-consistency relation p = ^, m,/L requires that x/,,,/, 
zkmax/E should be larger than 1 . So we get m/,,,^ : 



1/8 

^hub 



orz 



Enij^^f^/kmax- Then, the effective fugacity for other nodes is 
written as x, = h/kc where kc is the crossover degree defined 
as 



kr = 



'^max 



'hub 



(7) 



The effective fugacity ranges fromxm,„ = (ko/kc). 

The self-consistency relation p = Y.i mt/L is used to find the 
value of kc- It takes a different form depending on whether 
x„„„ is larger than 1 or not. First, let us assume that x,„,„ > 1, 
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TABLE I: Scaling of the crossover degree scale kc, the mean occu- 
pation number, and m/,„/,. 
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kc 
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0(L/lnL) 


8<5c 
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or equivalently, ko > . Then, m, ~ x/ for all nodes, and 
the self-consistency relation becomes as p = k^l^ /k^^ . It is 
compatible with the assumption ko > kc only when fe^/^ is non- 
divergent. Note that = J dkk^/^PdegXk) - J^"' dkk^/^-y. 
It is non-divergent when 5 > 5c where 



In that case, the mean occupation number is given by m,- = 
ck^^ with a constant c = p/k^/^. The hub has the largest mean 
occupation number, but it scales sub linearly with L as mi„ib ^ 
kmax ^ L^^/^ with 8c/8 < 1. Therefore we conclude that the 
system does not show the condensation for 5 > 5^. 

Let us assume the opposite situation where jCmm <C 1 or 
ko ^ kc- Then, the self-consistency relation is decomposed 
into two parts as p = p„ + ps for nodes with ki < kc and 

ki > kc, respectively. According to Eq. (|6j, they are given by 

1 /6 

Pn = Lki<k,Xi/L and p, = Lki>k,x/ /L. The summands have 
no singularity, and one can write in the °o limit as 

Pn = dkkPdegXk) 

Jko 

Ps - k;'^^ dkk'i^Pd,g,{k) . 

Jkc 

Note that /^^'' dkkP^eg, (k) is smaller than the mean degree k, 
hence is finite. So p„ vanishes as ^ k^^, which yields that 
Pj = p and kl^^ = j^^""'dkk^l^PA^^{k)/p - J^'""' dkk^/^-f. It 
is compatible with the assumption that kc ^ ko only when the 
integral is divergent, which requires that 5 = 5c or 8 < 5c. In 

either case, we obtain that kc ^ [lnA:,„„x]^'' or kc ^ k]nax^^\ re- 
spectively. The occupation number distribution is then easily 
found using Eqs. (|6j and (Q, which is summarized in Table|I] 
For 5 < 5c, the mean occupation number, as a function of the 
degree, satisfies the scaling form 

mk - G^{k/kc) , (9) 

where the scaling function behaves as Gg{y ^ 1)^3' and 

The maximum degree scales as k,„ax ^ L^'^, which gives 
kc ~ L^'^^ for 5 < 5c. Using this in TablelUfor 5 < 5c, we find 
that the nodes with k > kc are occupied by macroscopic num- 
bers of particles as mkykc ^ I^{k/k,„axy^^ = Further- 
more, the whole fraction of particles is distributed on those 




FIG. 1: Scaling of the mean occupation number at 5 = 0.0 in (a), 
0.2 in (b), 0.5 in (c), and 0.8 in (d). Data from different network 
sizes are distinguished with different symbols (o forL = 1000, □ for 
L = 2000, o for L = 4000, and A f or L = 8000). The slope of solid 
lines is 1 and that of dotted lines is 1 /5. 

nodes, which is clear from p„ ~ A:^' ^ L^'^'^^'. Therefore, 
we conclude that the system displays the complete condensa- 
tion for 5 < 5c. At the marginal case 5 = 5c, the complete 
condensation also occurs with a logarithmic correction (see 
TableU). Note that the case with 5 = corresponds to that with 
kc = kmax SO that whole fraction of particles are condensed on 
to the hub. The results at 5 = are also summarized in TableU 
The analytic results are confirmed by Monte Carlo simula- 
tions on the archetypical Barabasi-Albert (BA) network 1 1^ 
with Y= 3 and (k) — 4. It has 5c = 1 /2. The network sizes are 
L = 1000,2000,4000,8000 and the total density of particles 
is fixed to p = 2. Figure ^ shows the stationary state occu- 
pation number of each node as a function of the degree. For 
5 < 5c, the data follow the scaling form in Eq. (|9}. It is clearly 
seen that only the hub is the condensed node at 5 = 0. At 
5 — 0.8 > 5c, the condensation does not occurs and m ^ k'^l^ 
for all nodes. 

In the context of transport on a SF network, one can avoid 
the complete condensation if the transport capacity is de- 
signed to be greater than p{n) = n^' . The threshold value for 
5 is determined by the network structure, i.e., 5c = 1/ (y— 1). 
When the complete condensation occurs, the high degree 
nodes are flooded by macroscopic numbers of particles. The 
condensation may slow down the dynamics leading to conges- 
tion. We will study the dynamic aspect of the condensation. 

Although the ZRP is solvable as far as the stationary state 
properties are concerned, little has been known for its dynam- 
ical properties 1,10. 13.1 . We investigate relaxational dynamics 
of the ZRP on the BA network numerically. Initially particles 
are distributed uniformly over the network. Observed is the 
time evolution toward the stationary state. 

We propose that the quantity I{t) = (n,(f) /N)^ be use- 
ful in studying the condensation dynamics. It will be called 
the inverse participation ratio (IPR). Roughly speaking, the 
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FIG. 2: Relaxation dynamics at 5 = 0.2 in the BA network, (a) Scal- 
ing plot of ///stat. versus t/L~ with z = 0.83. (b) Scaling plot of 
mi./f'/('"^' versus /t/?^'/''"^' in the BA network of L = 10^ nodes. 



IPR is given by the inverse of the number of nodes whose 
occupancy is non-neghgible. The IPR in the stationary state, 
^stat. = lim,^oo/(f), can serve as the order parameter for the 
condensation, for it is nonzero only in the condensed states. 

We performed the Monte Carlo simulations on the BA net- 
works. In the condensed phase at 5 < 5^ = 1 /2, we found that 
the IPR converges to a finite stationary state value /stat. and 
that the approach to /jtat. becomes slower as L increases. It 
turns out that the relaxation time follows the power law 



(10) 



As an evidence for the argument, we measured the mean 
occupation number in the transient times f <C X in the BA 
network. We find that at each time step t there exists a degree 
scale k — k{t) such that all nodes with k <k behave as in the 
stationary state in the smaller network with the maximum de- 
gree k. Namely, mi^^^f-^ at time f is given by the formula in 

Table m with k„u,x and L being replaced with k and L' ~ k^^^^, 
respectively. Then, the scaling argument can be validated by 
showing that t - L'(i-^) - ^(i-S)/8^, i.e., k{t) - f8./(i-5). We 
confirm the scaling numerically with the scaling plot of nik/L' 
versus k/k at various times in Fig.|2](b), where all data col- 
lapse at A: ~ ^. We have also investigated the relaxation dy- 
namics through the IPR and m^- in various SF networks with 
different values of y and 5, and found that all numerical results 
are consistent with z = 1 — 5 |^| . 

In summary, we have studied the condensation phenomena 
in the ZRP with the jumping rate p{n) — on SF networks 
with the degree exponent y. We have shown analytically that 
the complete condensation occurs when 8 < 8^ = l/(y — 1). 
We have also shown numerically that the condensed phase 
with 8 < 8c has the slow relaxation dynamics, which is char- 
acterized by the power-law scaling of the the relaxation time. 

Acknowledgment: This work was supported by the KOSEF 
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with a dynamic exponent z- Measuring the value of i from 
I{t —i)— fl/stat. with a constant a — 0.8, we obtain that z — 
0.50, 0.66, and 0.83 for 8 = 0.5, 0.4, and 0.2, respectively. 
Figure|2l(a) confirms the power-law scaling at 8 = 0.2. On the 
other hands, the relaxation time does not grow algebraically in 
the non-condensed phase 1 15]. Therefore we conclude that the 
dynamics becomes slow when the condensation occurs. 

We present a scaling argument for the power law T ~ LF. In 
the periodic lattices the relaxation time is dominated by the 
process where macroscopic condensates at two sites merge 
into a single one |10]. Similarly, we assume that it is domi- 
nated by the particle transfer process between two nodes with 
the largest occupancy, that is, the hub and the second largest 
hub. Close to the stationary state, they are occupied by (L) 
particles and a finite fraction of them should be transfered to 
the hub. At each time step, (L^) particles jump out of the 
nodes and drift toward the hub. We further assume that the 
drift time does not scale algebraically with L, which might be 
justified by the small world property of the SF networks f5"l. 
Under these assumptions, we obtain that the transfer process 
takes 0{L/L^) time steps, hence i ~ U with z = 1 — 8. Note 
that it is consistent with the numerical results obtained from 
the IPR up to 10%. 
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